Interacting dark matter contribution to the Galactic 511 keV gamma ray emission: 
constraining the morphology with INTEGRAL/SPI observations 



Aaron C. Vincent*, 1 Pierrick Martin^, 2 and James M. Cline^ 1 

department of Physics, McGill University, 3600 Rue University, Montreal, Quebec, Canada H3A 2T8 
z Institut de Planetologie et d'Astrophysique de Grenoble, BP 53, 38041 Grenoble cedex 9, France 

We compare the full-sky morphology of the 511 keV gamma ray excess measured by the IN- 
TEGRAL/SPI experiment to predictions of models based on dark matter (DM) scatterings that 
produce low-energy positrons: either MeV-scale DM that annihilates directly into e + e~ pairs, or 
heavy DM that inelastically scatters into an excited state (XDM) followed by decay into e + e~ and 
the ground state. By direct comparison to the data, we find that such explanations are consistent 
with dark matter halo profiles predicted by numerical many-body simulations for a Milky Way-like 
galaxy. Our results favor an Einasto profile over the cuspier NFW distribution and exclude decaying 
dark matter scenarios whose predicted spatial distribution is too broad. We obtain a good fit to the 
shape of the signal using six fewer degrees of freedom than previous empirical fits to the 511 keV 
data. We find that the ratio of flux at Earth from the galactic bulge to that of the disk is between 
1.9 and 2.4, taking into account that 73% of the disk contribution may be attributed to the beta 
decay of radioactive 26 A1. 



1. INTRODUCTION 

The 511 keV gamma ray line observed by the INTE- 
GRAL/SPI experiment is consistent with the annihila- 
tion of ~ (1.5 ± 0.1) x 10 43 low-energy positrons per sec- 
ond in a region within ~ 1 kpc of the galactic center 



(GC), in addition to a fainter ((0.3 ± 0.2) x 10 



43 e + 



- 1 ) 



disk-like component that extends along the galactic plane 
[l[. The line is mostly due to parapositronium annihila- 
tion of thermal or near-thermal positrons @, II . The 
absence of 7 rays from e + annihilations in flight implies 
that the positrons are injected with energies less than 
~ 3 MeV [4J ■ No astrophysical source has been proven to 
yield such positrons with the required concentrated and 
approximately axially symmetric spatial distribution. 

Among conventional sources, radioactive ejecta from 
stars, supcrnovae and gamma-ray bursts can produce a 
large enough rate of positrons through (3 + decay, but 
their spatial distribution is not sufficiently confined to- 
ward the bulge: they predict a ratio of bulge to disk 
luminosities B/D < 0.5, whereas observations demand 
B/D > 1.4. Other proposed mechanisms also suffer from 
this problem. In addition, positrons from pair creation 
near pulsars or from p-p collisions associated with cosmic 
rays or the supermassive black hole tend to be too ener- 
getic. Low-mass X-ray binaries have received attention 
as a possible source, but these also do not give rise to 
large enough B/D [5| . A comprehensive review of these 
sources and the challenges they face is given in 

Dark matter (DM) interactions have the potential to 
explain the observed excess, either through direct annihi- 
lations of light (~ few MeV) DM particles into e + e~ pairs 
, or by the excited dark matter (XDM) mechanism, in 
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which excited states of heavy DM (%) are produced in 
X-X collisions, with subsequent decay of the excited state 
into the ground state and an e + e~ pair 0, The lat- 
ter scenario has the theoretical advantage that the DM 
mass is relatively unconstrained, requiring only that the 
splitting between the ground and excited states be less 
than a few MeV. 

XDM as an explanation for the INTEGRAL/SPI 511 
keV excess came under greater scrutiny in recent years af- 
ter it was proposed [Io| that nonabelian DM models could 
naturally have small ~ MeV mass splittings and simul- 
taneously explain additional recent cosmic ray anomalies 

11 , [l2| as well as hints of direct DM detection [l3[ . Ref. 

14 found that it is not possible to get a large enough 
rate of positrons for 511 keV emission in the nonabelian 
models that require production of two e + e~ pairs (one 
at each interaction vertex). However, the original model 
of @ can give a large enough rate [lj| since only one 
such pair need be produced, which is energetically eas- 
ier. Moreover, variant models involving metastable DM 
that scatters through a smaller mass gap 0, [l7| also 
give a large enough rate, and are largely free of threshold 
velocity issues. 

The aforementioned studies focused primarily on 
matching the overall rate of positron production, either 
ignoring morpholo gica l constraints or estimating them in 
a rough way. Ref. [18| is the only rigorous analysis with 
respect to dark matter models, done at a time when rel- 
atively little data had yet been accumulated. More re- 
cently, ref. [l^| carried out a study of DM predictions for 
the 511 keV angular profile, but comparing to a previous 
fit to the observed shape [2(j] rather than directly to the 
data. 

Our purpose in the present work is to improve upon 
these earlier papers by testing the DM model shape pre- 
dictions directly against the most recent INTEGRAL 
data. We will then examine how these DM models com- 
pare to the phenomenological models obtained in previ- 
ous studies, such as 12ll [23| , where the 511 keV celestial 
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signal is represented by analytical shape functions with 
several free parameters. As we will see, an interesting 
feature of the DM models is that their predictions de- 
pend on far fewer parameters and they can thus be a 
more attractive candidate if they are shown to provide 
as good a fit as the phenomcnological parametrizations. 

In the remainder of the paper, we first present the 
known sources of positrons in the galaxy, before dis- 
cussing our procedure for modeling the 511 keV sky in 
Sections [3] and |U We give our main results, along with 
the details of our fitting procedure, in Section [5] and 
briefly discuss the implications of this study in Section [5] 



2. KNOWN BACKGROUNDS 

In order to correctly model the possible contribution 
to the 511 keV signal from DM scattering, it is necessary 
to subtract from the data the contributions from known 
sources of low-energy positrons. They can be produced 
from j3 + decay of 26 A1 expelled from massive stars, as well 
as from 44 Ti and 56 Ni produced in supernovae. These 
contributions should be correlated with the stars in the 
galaxy, thus contributing dominantly to the disk compo- 
nent of the observed signal. 

The contribution of 26 A1 can be more directly assessed 
than that of the other radio-isotopes. During 26 A1 decay, 
the de-excitation of the resulting 26 Mg nucleus produces 
a gamma ray signal at an energy of 1809 keV whose mag- 
nitude and morphology has also been mapped by INTE- 
GRAL/SPI (H)]. Since each decay produces a positron 
and an 1809 keV photon, one can unambiguously deter- 
mine the fraction of the 511 keV signal originating from 
26 A1. Ref. [J showed that it accounts for roughly half 
of the disk component of the 511 keV signal, and we will 
confirm this. The contribution of 44 Ti and 56 Ni positrons 
cannot be evaluated in that way because of their shorter 
lifetimes. A corollary is that positron escape from su- 
pernova and their remnants can be a serious issue, and 
prevent the determination of positron injection rate di- 
rectly from the isotopes yields [24| [H| . Estimates of the 
isotopes production in stars and of positron escape frac- 
tions suggest that it should make up most of the remain- 
ing disk emissivity 0, Q . 



3. DARK MATTER HALO PROFILE 

Many-body simulations of the formation of galactic ha- 
los by collapsing dark matter particles predict a triaxial 
halo (see for example [26[), which however becomes more 
approximately spherical near the galactic center when 
the effects of baryons are taken into account [27j |. For 
simplicity we will consider the halos to be spherically 
symmetric in most of the present work, although we will 
show that adding a realistic degree of oblateness does not 
significantly alter the fit. To further constrain the shape 
of the halo we will refer to results of the Via Lactea II 



simulation |28j, which modeled the collapse of a Milky 
Way-sized (2 x 10 12 M ) collection of over 10 9 particles. 
We chose Via Lactea II because it was specifically geared 
towards the study of the dark matter halo of the Milky 
Way. Among the many known parametrizations of the 
radial mass-energy density distribution, two have been 
especially successful at parametrizing results of recent 
simulations. These are the Einasto profile 



p(r) = p s exp 



- 1 



(1) 



and the generalized Navarro-Frenk- White (NFW) profile, 



P(r) = Ps 



2 3- 7 



(r/r s )7(l + r/r s ) 3 -T' 



(2) 



In both cases r is the galactocentric radius, while r Sl a 
and 7 are parameters fit to N-body simulation results. 
The main galactic halo of the Via Lactea II simulation 
may be fit to an Einasto profile with r s = 25.7 kpc and 
a = 0.17, or to an NFW profile with r s = 26.2 kpc 
and a central slope of 7 = 1.2 [29[. The overall density 
normalization p s can be computed from the local dark 
matter density which we take to be pq = 0.4 GeV cm~ 3 
[30l ] at the sun's position r = 8.5 kpc f3l| . 



4. DM AND THE 511 KEV SKY DISTRIBUTION 




was already 
191. for com- 



Although the decaying DM scenario [3 
shown to be highly disfavored in refs. [1 
pleteness we will retest it in the present work. The flux 
of 511 keV photons from an e + produced in the decay of 
a metastable DM particle \ of mass m x is 



d$ = 2(l-0.75/ p )^ / 

47T h.o.. 



(3) 



The integral is along the observer's line of sight 
parametrized by £, r is the lifetime, p(t) is its position- 
dependent density and f p = 0.967±0.022 is the positron- 
ium fraction Q . It corresponds to the global probability 
that a given e + e~ annihilation take place via positronium 
formation. The latter can occur in the triplet state ortho- 
positronium (o-Ps) or the singlet state para-positronium 
(p-Ps). To conserve angular momentum, only p-Ps may 
decay into two 511 keV photons. 

If the positrons are instead produced in a scattering 
or annihilation event, the observed flux takes a similar 
form: 



d<S> 



2(1-0.75/,)^ 



l.o.s 



1)17 



d£ (4) 



where (av) is the thermally averaged cross-section for 
annihilations or excitations of the DM particles that pro- 
duce e + e~ pairs. Henceforth we will use "scattering" 
as shorthand for either XDM scattering or annihilating 
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light DM, since both processes will look like ((4]) to an ob- 
server. The density-squared dependence of this integral 
means that the observed flux is much more concentrated 
in the galactic center than in the decay case; this is why 
scattering gives a much better fit to the observed shape 
than do decays. 

The forms (|3I4[) are only strictly correct if positrons 
annihilate close to where they were formed. Despite re- 
cent studies (33l |3~H the problem of positron transport 
in the interstellar medium cannot be considered as fully 
settled. In the absence of strong theoretical and obser- 
vational constraints, we will assume that positron trans- 
port is a small effect in the present investigation. We will 
briefly return to this issue in Section El 

Moreover, we have for simplicity assumed that (av) in 
(|4"| is independent of r, but this is not a good approx- 
imation for all models. In particular, for the standard 
XDM scenarios with a total energy gap SE > between 
the ground state and excited state(s), there is a threshold 
value for the relative velocity, Vt = 2 s/SE /m x , which ap- 



pears in the excitation cross section as av <~ <jq \A>' 
Q. Because the DM velocity dispersion vo(r) depends 
strongly upon r near the galactic center, this factor can 
then introduce significant r dependence into the phase- 
space average (crv). There are several situations where 
this is not important: MeV DM undergoing pure annihi- 
lations @, I35J. metastable XDM models where 5E <C m e 
or 5E < |l4l. [13] , and standard XDM models where 
m x > TeV, in which case Vt is small compare to vq(t). 
For XDM models with m x < TeV, a more detailed study 
should be done. 

In addition to the dark matter source of positrons, we 
included a disk component that models (3 + emission from 
radioactive isotopes including 26 A1 and 44 Ti, whose flux 
at earth is analogous to eq. ©; the combination p/(m x T) 
becomes a density per unit time h of positron-producing 
radioactive decays. We considered two density models 
for this component. The first is a Robin young stellar 
disk (YD) model P,H|, 



h YD (x,y, z) = h Q 



with 



1 2 2/2 

+ y +z je 



(5) 



(6) 



The fixed disk scale radius is Rq = 5 kpc and the fixed 
inner disk trucation radius is Ri = 3 kpc. We varied 
the vertical height scale zq = e/Ro between 50 pc and 
140 pc. (Ref. [i| used the 1809 kcV line to fit the 26 A1 
distribution to a YD distribution with zo = 125 pc.) For 
comparison we also took an old disk (OD) model: 



h OD (x,y,z) = ho 



0.25+^ 



0.25-H 



1/2-, 



(7) 

with Rq = 2.53 kpc, Ri = 1.32 kpc and a vertical height 
scale zq which was varied from 150 to 250 pc. 



5. RESULTS 

We tested our DM scenario against the INTE- 
GRAL/SPI data by a model-fitting procedure applied 
to about 8 years of data collected in an energy bin of 5 
kcV width centered around 511 keV. For this, a model 
for the sky emission is convolved by the instrument re- 
sponse function and fitted to the data simultaneously to 
a model for the instrumental background noise in the Ge 
detectors. 

Our fitting procedure is the same as the one described 
in section 4.2.1 of [l|. The likelihood L of a model as- 
suming a Poisson distribution of events in each of the N 
data bins is 



L 



N 

n 

i=l 



A"* 



(8) 



rii is the number of events recorded in bin i by the SPI 
experiment, and A^ = J2k a k s i + bi(f3) is the predicted 
number of counts per bin, including the background bi 
and the source = Yljfj-^i- The factor Rj is the 
instrument response matrix and fj is the intensity com- 
puted with the line-of-sight integrals. In our case, the 
sum over k has two terms: the dark matter term and the 
disk component. The coefficients and /3 are the scaling 
factors that are adjusted by the fit. The result of fixing 
the normalization udm is to fix (m x r x ) in the case of 
decay and (av) x m~ 2 for dark matter scattering. We use 
the maximum likelihood ratio test to estimate detection 
significances and errors. We calculate the log-likelihood 
ratio 



MLR = -2(lnX -InLi), 



(9) 



where L\ is the maximized likelihood of the model being 
tested, and Lq is the maximum likelihood of the back- 
ground model only, i.e., = 0. 

We compare the results of our DM models to the phe - 
nomenological description by Wcidenspointncr et al. [2l| , 
where the authors fitted two spheroidal Gaussians and a 
young stellar disk to the then-available four-year data 
set 1 . We have updated their analysis, using the cur- 
rently available eight-year data set and find an MLR of 
2693. Although non-nested models cannot be directly 
compared through the MLR, this serves as a the figure 
of merit for a model such as the dark matter ones to 
match, if it is to provide a competitive fit relative to the 
phenomcnological shape models. 

We performed two analyses, firstly fixing a and r s to 
values favored by Via Lactea II, using the yo ung disk 
model parameters favored by the 26 A1 analysis of [23[ , and 



The 8 degrees of freedom in the reference model are: the width 
and normalization of each Gaussian, the inner and outer disk 
truncation, the disk scale height and the disk normalization. 
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Table I: Summary of best fits to the INTEGRAL/SP1 data, with parameters fixed to results of the Via Lactea II simulation. 
This corresponds to r s = 26 kpc and a = 0.17 for an Einasto profile {T} or 7 = 1.2 for an NFW profile (|2jl. The disk component 
is the young disk (0) with zo = 125 pc. All-sky fluxes are in units of 10 -4 ph cm _2 s _1 , the lifetimes r are in seconds, and 
cross-sections (av) have units of cm 3 s _1 . We have highlighted the best fit scenarios in bold. 



Channel 


Profile 


MLR 


Disk flux 


DM flux 


DM lifetime or cross-section 


decay 


Einasto only 


2f39 




174.5 ± 3.5 


r x = 1.1 x 10 26 (GeV/m x ) 




Einasto + Disk 


2f94 


10.60 ± 1.42 


148.6 ± 5.1 


r x = 1.3 x 10 26 (GeV/m x ) 




Einasto only 


2611 




24.02 ± 0.47 


{av) x = 5.8 x 10- 25 (m x /GeV) 2 


scattering 


Einasto + Disk 


2668 


9.98 ± 1.32 


21.16 ± 0.59 


(erv) x = 5.1 X 10~ 25 (m x /GeV) 2 


Einasto (oblate) + Disk 


2669 


8.74 ± 1.31 


21.06 ± 0.61 


(erv) x = 4.9 X 10 -25 (m x /GeV) 2 




NFW only 


1602 




6.72 ± 0.17 


{av) x = 8.2 x 10~ 26 (m x /GeV) 2 




NFW + Disk 


2155 


26.45 ± 1.25 


4.90 ± 0.18 


(av) x = 6.1 x 10^ 26 (m x /GeV) 2 




Figure 1: Intensity skymap predicted by Einasto + disk 
model. The bulge component is due to emission from scat- 
tering or annihilating dark matter in an Einasto profile, and 
the disk component can be attributed to decay of radioactive 
species including mainly 26 A1. 



finding the overall normalizations of the disk and Einasto 
components that best fit the INTEGRAL /SPI data. As 
a second analysis, we varied the parameters a and r s of 
the Einasto profile, as well as the height scales z$ for both 
young and old disk populations. As we will show, adding 
these three extra degrees of freedom does not signifcantly 
improve the likelihood of the model, suggesting that the 
Via Lactea II parameters are a good fit for the scattering 
XDM or annihilating DM hypothesis. 

Table [J summarizes our main results. The dark matter 
halo parameters were set to those favored by Via Lactea 
II, for an Einasto (NFW) profile with r s = 26 kpc and a 
= 0.17 (7 = 1.2). We used the young disk model ([5]) of 
[23l ] . with the fixed scale height zo = 125 pc corresponding 
to the 26 A1 distribution inferred from 1809 keV line data. 
We considered both decaying ([3]) and scattering (0| dark 



Figure 2: Longitudinal dark matter profiles for the three dark 
matter models considered, including the disk component from 
radioactive isotopes. Fluxes are integrated over galactic lati- 
tudes — 15° < b < 15°. "Scattering" refers to either scatter- 
ing multistate dark matter or annihilating light dark matter. 
The solid magenta line is left-right averaged, reconstructed 
SPI data from [f|, taken from the skymaps of [37| . 

matter. The scattering scenario provided a consistently 
better fit (AMLR> 400), and the fit to the Einasto pro- 
file was significantly better than to the NFW profile (A 
MLR= 513). Motivated by the triaxial halo shapes men- 
tioned above [27|, we also examined an oblate Einasto 
profile with a semi-major axis ratio c/a = 0.8. This is 
denoted "Einasto (oblate) + disk" in Table |U While this 
reduced the required flux from the disk component, it did 
not produce any significant change in MLR. 

The best-fit lifetimes (cross-sections) of the XDM 
model in the decaying (scattering) scenario are presented 
in the final column of Table |U Figure [T] shows the all- 
sky map of the Einasto + disk best fit to the INTE- 
GRAL/SPI data, and Figure [5] shows the longitudinal 
profile of the three dark matter models (including disk 
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'5 10 15 20 25 

r, (kpc) 

Figure 3: Maximum log-likclihood ratio (MLR) obtained in 
the decaying dark matter + young disk scenario as a function 
of the Einasto halo parameters. The values favored by the 
Via Lactea II N-body simulation, labeled VL2, do not give a 
good fit to the INTEGRAL/SPI data and are far away from 
the favored region. 



parameters. The reduced x 2 of our dark matter model 
computed on a pointing basis is as good as that of the 
phenomenological model, with a value of 1.007. 

Letting r s , a and zq vary freely yields some improve- 
ment. Figure[3]shows a contour plot of the MLR obtained 
from the decay scenario J3]). The favored region in the 
lower-left corner, with an MLR of 2558, corresponds to 
an extremely cuspy DM halo that is quite far removed 
from realistic DM halo models. 

The equivalent picture for scattering DM is illustrated 
in Figure The overall best fit was found to be for 
a profile with a = 0.2, r s = 12 kpc and zq = 140 pc, 
with an MLR of 2673. However, this difference is only 
marginally significant. Indeed, by adding three degrees of 
freedom, such an improvement should happen by chance 
17% of the time due to statistical fluctuations in the data. 
We found that the young disk (YD) model consistently 
gave a better fit than the old disk (OD) model, and that 
adjusting zo over a range from 70 to 200 pc did not pro- 
duce any significant improvement in the MLR. Finally, 
we checked that choosing a closer value for the galacto- 
centric distance of R. — 8.2 kpc, as suggested by recent 
studes such as [3(| produced a negligible change in the 
fit (AMLR < 1). 




r s (kpc) 

Figure 4: Same as Figure [3] but with scattering dark matter 
(Q. The MLR obtained with the Via Lactea II parameters 
(white dot) is within AMLR = 5 of the best fit (r 3 = 12 kpc, 
a — 0.2), which means that the VL2 parameters likely cor- 
respond to the correct model if the scattering or annihilating 
dark matter hypothesis is true. 

components) in comparison with a reconstruction of the 
SPI data. This clearly illustrates how decaying dark mat- 
ter produces a profile that is far too flat, while an NFW 
distribution results in an unrealistic sharp central peak. 
Decaying dark matter in an NFW profile (not illustrated) 
displays a combination of these flaws. On the other hand, 
the scattering model produces MLR = 2668, which is not 
far below that of the best-fit phenomenological model, 
the latter having MLR = 2693 and six additional fitting 



6. DISCUSSION AND CONCLUSION 

We have made the first direct comparison of dark mat- 
ter predictions for the observed 511 keV spatial inten- 
sity distribution since the earliest data release of IN- 
TEGRAL/SPI. Our favored fit corresponds to a scat- 
tering excited DM or annihilating light DM model in an 
Einasto density distribtion ([T} with parameters fixed to 
the Via Lactea II results. We confirm previous anal- 
yses showing that decaying dark matter is ruled out 
due to its too-broad spatial distribution. After cor- 
rect normalization of the intensity, our best-fit model 
requires a cross section for \X to produce positrons of 
(av) x = 5.1 x 10~ 25 (m x /GeV) 2 cm 3 s _1 . If m x is in the 
10-1000 GcV range as favored by most WIMP models, 
this means (an) is in the interval [l0~ 23 , 10~ 19 ] cm s . 
The fact that this is far above the annihilation cross sec- 
tion of 3 x 10~ 26 cm 3 s _1 needed to get the observed relic 
density is not problematic, because the physical process 
required in these models is inelastic scattering to an ex- 
cited state rather than annihilation. 

Because we neglected r-dependence in the averaged 
cross section (crv), these results apply to upscattering 
XDM with hi gh masses m x > a few TeV, metastable 
XDM models [H, [l7[, and direct annihilation of McV 
DM. To cover the case of lighter XDM models, a more 
detailed analysis taking account of the radial dependence 
of the DM velocity dispersion in the Galaxy would be 
needed. We hope to return to this in future work. 

For light ~ MeV DM annihilating directly into e + e~, 
our required cross section is (av) ~ 10~ 31 cm 3 s _1 , which 
is too small to give the right relic density. This need 
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not be a problem; it only requires there to be additional 
stronger annihilation channels into invisible particles, for 
example dark gauge bosons [35| or dark neutrinos [38l |. 

There are two unknowns that could change our anal- 
ysis in significant ways. One is the distance by which 
positrons propagate between creation and annihilation. 
If it is larger than ~ 100 pc, it could alter the overall 
breadth of the spatial extent of the signal, as well as intro- 
duce deviations from axial symmetry, depending on the 
conditions of the interstellar medium in the bulge. Fur- 
ther observational evidence constraining the structure of 
ma gne tic fields (for example synchrotron emission stud- 
ies [33| ) will be needed to reduce these uncertainties. A 
second unknown is the degree of departure of the DM 
halo from spherical symmetry, which definitely occurs in 
iV-body simulations [27j. We showed that adding some 
oblatcncss had little effect on the fits, though the nature 
and extent of triaxiality near the galactic center depends 
heavily upon the inclusion of baryons in the simulations, 
a challenging field which is still in its early stages. We 
look forward to improvements in these studies that will 
help to constrain the theoretically expected extent of tri- 
axiality in the DM halo. 

We have confirmed the findings of previous studies con- 
cerning the disk emission. Given a young disk model for 
the distribution of 26 A1, the observed flux of 1809 keV 
gamma rays [23j translates into an expected 511 keV flux 
of (7.33 ± 0.89) xlO -4 ph cm _2 s -1 . This alone accounts 
for 73% of the disk component favored by our model. If 
similar amounts of 44 Ti are present in the Galaxy, there 
is no need for an extra component to explain the disk 
component of the 511 keV signal. On the other hand, 
simulations show that in addition to the DM halo, there 
may also be a DM disk. This would give an extra DM 
contribution to the disk component of the 511 keV emis- 
sion. However, there is as yet no direct evidence for a 



DM disk in our own galaxy [|(| IH| . 

It is worth emphasizing that only two degrees of free- 
dom were required to obtain the MLR of 2668 in the DM 
scattering/ annihilation scenario. This is in contrast to 
the 8 d.o.f. necessary to obtain an MLR of 2693 with 
one best-fit phenomenological model. A further advan- 
tage of the DM model is that it is motivated by particle 
physics and cosmology, and it has a concrete, calculable 
production mechansim for the excess electron-positron 
pairs. Our results are independent of the details of the 
DM model, so long as the scattering events lead directly 
to an e + e~ pair. 

We find these results to be encouraging for the dark 
matter interpretation of the 511 keV excess, an anomaly 
that was first seen in 1972 by balloon-borne detectors 
[42^ . We hope that the experimental hard X-ray / soft 
gamma-ray astronomy community will be motivated to 
consider a higher-resolution instrument that would be 
sensitive to the 511 keV region of the spectrum in the 
future. Such observations would help to shed more light 
on this intriguing possibility, which could be the first 
evidence for nongravitational interactions of dark matter. 
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